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Abstract 

We develop a Feynman rule for energy-level diagrams emphasizing their connections to the 
double-sided Feynman diagrams and physical processes in the Liouville space. Thereby we com- 
pletely identify such diagrams and processes contributing to the two-dimensional response function 
in the Brownian oscillator model. We classify such diagrams or processes in quartet and numeri- 
cally present signals separately from each quartet of diagrams or Liouville-space processes. We find 
that the signal from each quartet is distinctly different from the others; we can identify each peaks 
in frequency domain with a certain quartet. This offers the basis for analyzing and assigning actual 
two-dimensional peaks and suggests the possibility of Liouville-space-path selective spectroscopy. 
As an application we demonstrate an example in which two familiar homogeneous mechanisms of 
relaxation are distinguished by existence or non-existence of certain peaks on the two-dimensional 
map; appearance or disappearance of certain peak is sensitive to the coupling mechanism. We also 
point out some confusion in the literature with regard to inclusion of relaxation effects. 
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I. INTRODUCTION 



The use of ultrashort laser pulse to probe the properties of molecules has been propelled 
by the rapid advances in laser measurement techniques. Recently, two-dimensional (2D) 
vibrational spectroscopy has been actively studied, where the spectral properties of multi- 
body correlation functions of polarizability (2D Raman spectroscopy) @, || |, ^ || |7], |8], |||, [H| 
|TT| , [T2| , [L3|, [TJJ or dipole moment (2D infrared spectroscopy) [16], [T7|, are measured. 
The 2D technique provides information about the inter- and intra-molecular interactions 
which cause energy relaxations. (ITJI 120, ETL E2L 



Theoretically, optical responses of molecular vibrational motions have been studied 



mainly by either an oscillator model |24] or energy level model. |25| The oscillator model 
utilizes molecular coordinates to describe molecular motions. This description is physically 
intuitive since optical observables (dipole moments or Raman polarizability) are also de- 
scribed by molecular coordinates; the effects of relaxation, which are caused by interactions 
of the coordinate with some other degrees of freedom, are rather easy to be included. As 
long as the potential is harmonic or nearly harmonic, signals can be calculated analytically. 



26, E7L 28, 29 



On the contrary, the energy-level model employs the energy eigen functions of a molecular 
motion but is physically equivalent to the oscillator model. Accordingly, laser interactions are 
described by transitions between the energy levels; the optical processes, including the time- 
ordering of laser pulses, are conveniently described by diagrams such as Albrecht diagrams, 
P2| or double-sided Feynman diagrams. JI[ Although the inclusion of relaxation processes 
from physical insight is less intuitive and is restricted to some special cases, this model has the 



advantage in identifying peak positions of optical signal in frequency domain. |33|, |34], |35, |36 



The anharmonicity of potential and nonlinear mode-mode coupling are also easily taken into 
account. Phase matching conditions, which chose a specific Liouville path contribution by 
the configuration of Laser beams, |l| is also easy to take into account. In the oscillator 
model or molecular dynamical simulations, the phase matching condition can be done only 
after calculating entire response functions. f3"7fl 

The rate of increase in the number of diagrams, however, with the increase of laser 
interactions is severer in the energy-level model; this becomes serous practical problem for 
multi-dimensional spectroscopy, where many laser interactions are included. For example, 



3 



more than 16 diagrams are involved in the lowest order in the third-order anhamonicity 
in fifth-order Raman while all the diagrams can be represented by a single field-theoretic 
diagram in the oscillator model, p? 



In this paper we try to bridge the two complementary models by transferring some results 
obtained in the oscillator-model to the energy-level language. Although we lose the simplic- 
ity (e.g. small number of diagrams) we gain in an insight into optical processes; we can 
assign each peaks in certain optical or Liouville-space processes. The resulting energy-level 
Feynman rule for the oscillator system allows inclusion of relaxation mechanism in an ad 
hoc way. As an application, we compare two system with different damping constants. This 
example reveals that existence of certain peaks in 2D spectroscopic map sensitively depends 
on the relaxation model. 



II. INTERACTION OF ENERGY LEVEL DIAGRAMS 

We consider a molecular vibrational motion described by a single molecular coordinate 
Q. In the energy-level representation, the Hamiltonian is expressed as 

H Q = htt fata + (1) 

where a and a) are the creation and annihilation operators and 



Q = \ — — (a + a ] ) 
^ V 2Mtt K > 



(2) 



for the system with the mass M. The energy level of this harmonic system is given by E n = 
HQ n with Q n = (n + 1/2)Q for which we introduce the frequency difference f2 mn = Q m — Q n . 
If the system interacts with the laser field E(t), it is governed by the full Hamiltonian, 

[ H + uE(t) (IR) 
H(t) = < (3) 
I H + aE(t) 2 (Raman) 

where /x is the dipole for infrared (IR) and a is the polarizability for Raman spectroscopy. 
Both operators can be expanded as 

x = x + xiQ + ^x 2 Q 2 + ^j x 3<5 3 H , (4) 
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We consider the response function, which is pertinent to the 2D second-order IR (for non 
isotropic media) or the 2D fifth-order Raman spectroscopy, 



R®(T U T 2 ) 

= o(h-t 2 )6(t 2 



x (* 3 ) , j% (h) 



where x(t) is the Heisenberg operator of x for the non- interacting Hamiltonian H and 
(O) = Tr [poO] with p = e _/3H °/Tr [e _/3H °] (when we include the effect of dissipation at the 
level of Hamiltonian, Hq includes the bath Hamiltonian and the system-bath interaction). 
The operator x stands for \x (IR) or a (Raman). Generalization to the combined IR and 
Raman cases such as ([[/x (t 3 ) , \i (t 2 )] , a (^i)]) @> P?| ^ wm also be treated below. 

(Ti,T 2 ) for the harmonic system can be expanded in terms of Q by Eq. (|]). The 
leading order is given as 



i? (2) m,T 2 ) 



2 2 



(i?! + # 2 + R 3 ) 



(5) 



where 



i?i = <[[Q 2 (7i + T 2 ) , Q (TO] , Q (0)] ) 

#2 = <[[Q (T 1 + T 2 ),g 2 (TO],g(o)]) 

i?3 = <o^m+^),g(T 1 )],Q 2 (o)]>. 



with 



^3 ~ h — T 2 
t 2 — t\ = Ti, 



(6) 



A. Raman spectroscopy 

For the moment, we concentrate on the Raman case, i.e. ([[a (t 3 ) , a (t 2 )] , a (ti)]). Some 
of processes in Eq. (||) are represented by the energy-level (Albrecht-like) diagrams in Fig. 
[I]. The differences from the original Albrecht diagram are mentioned at the end of this 
section. Before explaining diagrams, let us review possible transitions by operators Q and 
Q 2 ; Q can cause a one-quantum excitation or de-excitation while Q 2 can result in a two- 
quantum excitation or de-excitation in addition to a zero-quantum transition. For example, 
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from |0) — > Q 2 |0) ~ ( flt ) 2 + aflt j |0), we see that by the action of the operator Q 2 , the 
ground ket state |0) can be changed into |0) (zero-quantum transition) or |2) (two-quantum 
excitation). In the same way, (2| can be brought into (0| (two-quantum de-excitation) or 
<2|- 

In the diagrams, time runs from the left to the right. Each pair of arrows stands for a 
Raman excitation. The pair with a wavy arrow signifies the Raman induction decay (last 
interaction); the first interaction occurs at t±, the second at £2, and the last at £3. 

The full description of a quantum state at a certain time requires both the bra state (n\ 
and ket state |m); at any time the state is fully specified by the Liouville state \m) (n\. In 
the diagrams, the excitation or de-excitation of the bra state is expressed by a pair of fine 
arrows while that of the ket state by normal ones. For example, the first interaction at ti 
of (i) and (ii) is a two-quantum excitation of the ket state while that of (iii) and (iv) is of 
the bra state. 

In the Liouville space, the diagram (i) is interpreted as follows. The system is initially in 
the ground (Liouville) state |0) (0|. The first interaction causes a two-quantum excitation 
of the ket state; |0) (0| — > |2) (0| at t±. The second interaction causes a one-quantum de- 
excitation, |2) (0| — > |1) (0| at t 2 . The last shows a one-quantum de-excitation, |1) (0| — > 
|0) (0| at t 3 . As a whole, we denote this as 

|0)(0| r |2)(0| r |l)(0| r |0)(0| (7) 

ii £2 £3 

The diagram (ii)-(iv) are interpreted as follows: 

|0)(0| r |2)(0| r |l)(0| r |1)(1| (8) 

ii t 2 13 

|0)(0| r |0)(2| r |0)(l| r |0)(0| (9) 

tl 12 tz 

|0)(0| r |0)(2| r |0)(l| r |1)(1| (10) 

tl t-2 13 

Note here that a pair of fine arrows always correspond to the excitation or de-excitation of 
the bra state. 

We define the population state by \n) (n\, while the coherence state by \n) {m\ with 
n 7^ m. We notice that, after the last interaction, in all of the above 4 diagrams, the system 
is always in a population state (|0) (0| or |1) (1|). In summary, a diagram does not vanish 
only when the final state is a population state (Theorem 1). This corresponds to the trace 
operation in the definition of the response function. 
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FIG. 1: Energy- level diagrams of RW{T 1 ,T2) for Raman processes. 



In this paper, we simplify the original Albrecht diagrams for comparison with the 
Liouville paths. The main differences are the following: (1) we use always the same horizontal 
lines regardless of ket or bra states; it is not the case in the original Albrecht diagrams and 
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(iv) 







FIG. 2: An energy-level diagram of i?( 2 )(Ti, T2) for IR processes. 

(2) time runs always from left to right in our representation while the direction for the bra 
and ket states are the opposite in the original version. Our representation is somewhat 
simpler in that a single diagram in ours corresponding to several diagrams in the original 
version. 



B. IR and IR- Raman spectroscopy 

IR processes appearing in the IR response function, ([[// (t 3 ) , \i (t 2 )} , \x {ti)]), correspond- 
ing to Fig. |l|-(iv) is described in Fig. [| each quantum transition is represented not by a 
pair of arrow but an arrow. Note Raman and IR processes can be equivalent theoretically 
at this level of description, although even orders of IR processes, such as second-order IR 
signal vanish except in anisotopic media, such as adsorbed molecules on metallic surface. 



41j This situation can be overcome by mixing the IR and Raman processes. |38j By using 
narrow-band lasers (two IR excitation pulses followed one probe pulse which create Raman 
signal) Zhao and Wright demonstrated such experiment. 0, |(| 

As an IR-Raman spectroscopy, we consider the response function, 
([[// (£3) , /i (£2)] , ck (ii)]), for example. A diagram corresponding to Fig.^-(iv) is shown 
in Fig. [| Raman and IR transitions are represented by a pair of arrows and an arrow, 
respectively. Diagrams corresponding to the other IR-Raman response function such as 
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FIG. 3: An energy-level diagrams of R( 2 \Ti, T2) for IR-Raman processes, ([[// ^3) , \i fa)] , a fa)]). 
([[/i fa) , /j, fa)] , a fa)]) can be described in a similar manner. 

III. ENERGY-LEVEL DIAGRAM AND DOUBLE-SIDED DIAGRAM 

We can represent processes in the Liouville-space in a different way by the double-sided 
Feynman diagrams. The diagrams in Fig. |] are the translation of the diagrams in Fig. [I], [2|, 
or 0. In the double-sided diagrams, time runs from the left to the right (as in the energy-level 
diagram). The horizontal lines, however, are always two in number, the upper and lower 
line. The former represents the ket state while the latter the bra state. The single circle 
stands for a one-quantum transition, while the double circle for a two quantum transition. 
The quantum number of the bra and ket states is denoted explicitly in the diagram. 

It is noted that there are some differences of diagrammatic notation among articles. For 
example, in some literature, the quantum transition is not represented by circles but arrows. 
In other one, diagrams are rotated by 90 degree so that the time runs from the bottom to 
the top. 

In general, as seen below (VI. A), the double-sided diagram is convenient for enumeration 
of all possible diagrams while the energy-level diagram is for understanding the physical 
process. 
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FIG. 4: Double-sided Feynman diagrams of R( 2 \Ti, T2). 

IV. FEYNMAN RULES FOR THE DIAGRAMS 

We have introduced several way to represent optical processes as in Figs. |l]-§|. It is em- 
phasized here that the interpretation in terms of the Liouville-space state \m) (n\ is unique 
except for what Xk implies. Accordingly, we can develop a universal rule to write down 
analytical expressions from diagrams via the interpretations (such as Eqs. (|7|)- (|10D ) in the 
Liouville-space; the derivation is a straightforward exercise in elementary quantum mechan- 



ic) 



ics and would be discussed elsewhere. It can be summarized in the following way. We 
associate with each interaction (originating from the interaction Q k /k\) at a certain time or 
each propagation for a certain period one of the following factors: 



interaction (n > 0) 


factor 


\m) — > \m + n) 


Ix^ (m + n Q fe /k\ 


m) — > \m — n) 


(m — n\ Q \m) jk\ 


(m\ — * (m + n\ 


— \%k ( m \ Q k \ m + n) jk\ 


(m — > (m — n | 




remark 


omit ± | for the last interaction 



propagation (t > 0) 


factor 




m) (n 


for t 





By multiplying all the factors and putting another factor 1/2 to avoid double-counting 
(see Theorem 2 below), we obtain an analytical expression of the corresponding diagram 
(Feynman rule). Here, we have introduced ( mn and T mn (> 0) to describe relaxation; the 
difference of frequency modified due to the relaxation is defined by ( mn = (m — n)( while the 
relaxation constant T mn for the state \m) {n\ possesses the symmetric property, T nm = T mn , 
which is a necessary condition for a consistent theory (see below Eq. (|T2"D). Without 
dissipation, ( mn — > Q mn = (m — n)Q and T mn — > 0. In the Brownian oscillator model with 
the damping constant 7, the corrected frequency ( is given by ( = \Jvi 2 — (7/2) 2 . ^ 
The expression for T mn in this model shall be discussed below. 

By definition, the propagation period implies the time between two interactions. This 
excludes the periods from tj to ti and from t 3 to tp in diagrams in Figs. |l]-[| (or, say, in Eq. 
(0)"(HO)) because there is no interaction at tj or tp; we associate the unity for these special 
period. 

Let us apply our rule without relaxation (r m „ = 0, ( mn = Q mn ) to a diagram or a 
Liouville-space path. As the first example, we consider the diagram (i) (of Fig. [j] or 
We have only 2 separate propagation periods by definition. In the first period from ti to £2 
the system is in the state |2) (0| and thus we have the factor e - ^ 20 ^ 2- * 1 ) while for the last 



-ifl 10 (t 3 -t 2 ) . 



period from £2 to £3 the system is in the state |1) (0| and we have the factor e 
total we have the propagation factor, e~ lQ2oTl ■ e~ lQl ° T2 , where we have used the relation 
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In addition, as the result of the three interactions, we have other factors, ix 2 (2| Q 2 |0) /2 • 
j:Xi (1| Q |2) ■ xi(0|Q|l) = (i2l?r2* ri ) 2 X2 (Note nere the relations, Eq. (fj) as well as, 
a \n) = i/n \n — 1) , and a''" |n) = yn + l \n + 1)). In summary, the process in Eq. (|7]) or the 
diagram (i) is given (with the extra factor 1/2 associate with double-counting) by 



t%. Although the last interaction at £3 is that for the bra state (expressed by the fine arrows 
and different form (i)) the factors for this last interaction is the same with that of (i) by the 
above Feynman rule; there is no sign differences between bra and ket states (only) for the 
last interaction. In summary we have 



In general, we have the following theorem, which is related to the double counting: The 
diagrams different only by the side of the last interaction (bra or ket side) have the same 
contribution (Theorem2). 

The process in Eq. @ or in the diagram (iii) can be estimated in a similar manner by 
the above Feynman rule: 



Note here that the sign in front of i/h is minus because of the interactions on the bra state 
(fine arrows). From t\ to t 2 , the system is in the states |0) (2| and |2) (0| in (iii) and (i), 
respectively; these two states are the complex-conjugate of each other. From t 2 to t 3 , the 
state of (iii) (|0) (1|) is again in the complex-conjugate state of (i) (|1) (0|). Accordingly, (iii) 
given in the above is the complex conjugate of (i), i.e., (iii) = (i)*. Diagrammatically, in (iii) of 
Fig. m, all the normal arrows in (i) are replaced by the fine arrows. In general, The complex- 
conjugate diagram is obtained by interchanging all the normal and fine arrows (Theorem 3). 
In the double-sided Feynman diagrams, instead, The complex- conjugate diagram is obtained 
by interchanging the circles on the upper and lower lines (Theorem 3'). 

The diagram (iv) is the complex-conjugate diagram of (ii) because the fine and normal 
arrows are interchanged, i.e., (iv) = (ii)*. We can also verify the relation, (iii) = (iv), from 
the above Feynman rule with reconfirming Theorem 2. 




(11) 



The process in Eq. (||) or the diagram (ii) (of Figs. |T| or |j) is different from (i) only after 



11 = 1 
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V. TEMPERATURE EFFECT AND INITIAL STATE 



In the above, we have assumed the system is initially in the ground state |0) (0|, which 
is usually justified for high frequency vibration modes at a room temperature. For high 
temperatures or low frequency modes, however, excited states \n) (n\ are initially populated 
according to the Boltzmann factor. In general, we have to estimate all the possible processes 
assuming that the system is initially in the population state \n) (n\ using the above mentioned 
rule, and then summing up with respect to n with the Boltzmann factor e~^ En / e~@ En 
(in the case without dissipation); this completes our Feynman rule. 

Even if we take into account the contribution from general initial state \n) (n\, however, 
in the (fully-corrected) Ohmic Brownian oscillator model, we still have the same result with 
above as shown in the previous literature. This is the reflection of the relation 

(n\X\n) = (0\X\0) 

where X is some special combination of operators (This could be directly checked by la- 
borious calculation by using our Feynman rule). The fact that R^(Ti,T2) treated in this 
paper is independent of the temperature and thus we can obtain a finite temperature re- 
sult even if assuming that the system is initially in the ground state is by no means trivial 
but established by other calculation methods. |27| This implies, for example, that the de- 



pendence on n of the analytical expression corresponding to Fig. |5] cancels out with some 
other diagram. When the damping mechanism other than (fully-corrected) Ohmic Brown- 
ian oscillator model, our results presented below might be interpreted as an high frequency 
approximation, i.e. Ml^> kT. 



VI. LIOUVILLE-SPACE QUARTET 

The four diagrams (i)-(iv) in Figs. [I] and |] are a special set in the sense that we can 
obtain the other three, starting from one of the quartet. 

In the energy-level diagram, we obtain the second by changing the last interaction by 
using one of the following rule (depending on the last interaction of the starting diagram); 
(1) the ket excitation to a bra de-excitation, (2) the ket de-excitation to a bra excitation, (3) 
the bra excitation to a ket de-excitation and (4) the bra de-excitation to a ket excitation. The 
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FIG. 5: General process corresponding to Fig. [I] (i) 

remaining two diagrams are the complex-conjugate diagrams of the previous two diagrams 
where the conjugates are obtained by interchanging the fine and normal arrows. 

In the double-sided diagram the second diagram is obtained by lowering or raising the 
last circle. The remaining two is by interchanging lower and upper line with circles. 

As seen before, the corresponding analytical expressions of (i)-(iv) have the relations, 
(i) = (ii), (iii) = (iv), and (i) = (iii)*. The sum of the quartet is always real: 

(i) + (ii) + (iii) + (iv) = 4Re[(i)] = 4Re[(n)] 

where n — i, ii, iii or iv. Taking the real part of Eq. ([11]) we have an expression for the 
quartet, 

^ = --^^cos(2QT 1 + fiT 2 ) (12) 
Z {Mil) 

In terms of the interpretation in the Liouville space in Eqs. (]7j)-(|T(J|), all the processes 
posses a common property; the two-quantum coherence (|2) (0| or |0) (2|) is realized for T\ 
while the one-quantum coherence (|1) (0| or |0) (1|) for T 2 ; we denote this as: 

|2><0|-|1)<0| and |0> (2| - |0> (1| 

This is reflected by the factor cos (2QT\ + fiT 2 ) in Eq. (|T2"D . 

We notice that in the case with damping if T mn (and ( mn ) were not symmetric, Vj could 
not be real; the symmetric property of T mn is required for the response function to be real. 
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FIG. 6: All possible quartets for R ( - 2 \T 1 ,T 2 ). The square bracket implies that four diagrams 
are collectively represented. For example, the first diagram in the energy-level diagram for R(l) 
corresponds to not only (i) of Fig. [I] (which is explicitly written in the bracket) but other three 
diagrams (ii)-(iv) of Fig. |]. 

A. Quartets representation: all possible quartets for R <y2 \T\,T2) 

We show six quartets R(l)— R(6) in Fig. |]m the double-sided representation. The square 
brackets imply the quartet; only one of the quartet is explicitly written in the bracket. For 
example, R(l) of Fig. |5] collectively stands for (i)-(iv) of Fig. f|. 

In Fig. on the right side, ten quartets in the energy-level representation are given; 
some quartets in double-sided representation corresponds to not one but two quartets in the 
energy-level representation. For example, R(l) contains contribution I and I', while R(3) 
contains only A2. 
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These six quartets R(l)— R(6) in Fig. |6] exhaust all possible contribution to the right- 
hand side of Eq. (|5|); there are 3 ways for the position of the double quantum transition 
(double circle) and there 2 3 ways to put the three (including one double circle) circles upper 
or lower line, which leads to 3 ■ 2 3 double-sided Feynman diagrams in total. These 3 • 8 
diagrams can be divided into 6 quartets that have been shown. We understand here that 
the double-sided diagram is convenient for enumerating all possible diagrams. 



B. Estimation of quartets 



The analytical expression of quartet II is given via our Feynman rule: 



II = 4 Re 



2 / t \ 2 

i \ ( a \ 



h 



\2M() 



2 

x Xl%2 . e -* 2 C3i-r 20 Ti . e -i(T 2 -r 21 T 2 



(13) 
(14) 



where the analytical expression in the square bracket has been derived from the diagram 
explicitly drawn in the bracket in Fig. ^| (in the presence of dissipation). For example, the 
propagator e~* 2< » Tl-r2oTl and e~^ T2 ~ r2lT2 come from the propagation of |2) (0| and |2) (1|, 
respectively. 

In this way we obtain the expression: 



R {2) (7i, T 2 ) = I + II + A + B + C + D1+D2 



(15) 
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with 



I 

II 

A 
B 
C 
Dl 
D2 



x\x 2 



2 (M(Y 



x\x 2 



-r2o7i— rioT2 



cos (2CT! + CT 2 ] 



2 (MC) 

2 



2(MC) 



. e -r2oT!-r 21 T 2 cog + 
e -rioTi-rio2b cos (^ + ^ Ta ) 



2(MC) 
2(MC) 

„2 



1 X x X 2 c _ TioTi _t ooT2 



4 (MC)' 
3 x\x 2 



cos (CTi) 



A (M() 2 

As for the derivation of this we remark: (1) Quartets I' and IF cancel out because F- 



1 x i x 2 FnnTi — FinT?, 
4 (MC) 



?e ^uu^i ^u^cos((T 2 ) and IF 



i (mc) 2 6 r °° Tl r ° lT2 cos (C^2) (The numerical factor 



1/4 can be understood from the first two-quantum transition associated with (0| Q 2 |0) cc 
(0| aa^ |0) = 1). (2) The sum A2+A1 reduces to A (The numerical factor for A2 (or Al) can 
be estimated by noting the second two-quantum transition (1| Q 2 |1) cx (1| aa^ + a^a |1) = 3 
(or (0|Q 2 |0) a (0[ aa+ 10) = 1)). 

It is worth while observing the relationships between analytical expressions and the sym- 
bolic interpretations of the remaining quartets: 



A: 


|1) (Di- 


+ |1)(0| 


and |0) (1| - 


-|o)(i| 


B: 


li) (Di- 


+ |1)(2| 


and |0) (1| - 


-|2)(1| 


C: 


li) (Di- 


- |2) (0| 


and |0) (1| - 


-> ID) (2| 


D2: 


ll) (Di- 


-|1) (1| 


and |0) (1| - 


-|1)(1| 


Dl: 


ll) (o|- 


- ID) (0| 


and |0) (1| - 


- ID) (0| 



That is, we can associate the state \n) (m| with £ njm and T mn . 

In addition, if we fully include the temperature effect by our Feynman rule with tracking 
all the possible processes, we could obtain the result given in Appendix B of | j33fl . 
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VII. DAMPING MECHANISM 



We can confirm that the well-known result for the Ohmic Brownian oscillator (BO) model 
(Ohmic implies that the system-bath coupling is in the bilinear form) is reproduced from 
Eq. ( |I5D by setting 

r J 7 for |n > ^ fits 

r„ m = < (16) 

I \n — m\^f/2 for \n) (m\ (n ^ m) 

where \m\ represents the absolute value of m. Actually, in the Brownian result, I+II should 
be zero, which is true if T 2 i = r 10 , while D1+D2 should be — 2-D1, which is true if Tu = T 00 ; 
T mn in Eq. flTB] ) satisfies these requirements. 

The cancellation of I and II is one of the feature of the Brownian result. Another feature 
is that the state |0) (0| decays with the relaxation constant 7/2 which is the same as that 
for |1) (1|. These characteristics have intrigued some controversy as mentioned below. 

The relaxation constant for the same Ohmic model within the lower level approximation, 
i.e., at the level of the Fermi's golden rule with a somewhat ad hoc approximation (see 
below), given by |33], [44 



r m „ = — - — 7, (17) 

which is also simple but incompatible with the above two requirements (]?2i = Tio and 
Tu = Too)- With this relaxation constant, I and II survives, for example. (In addition, 
there is no frequency shift (£ mn — > VL mn ) in this finite-order approximation). 

The frequency shift and appearance of the absolute value (\n — m|), which is non- analytic, 
in the off diagonal relaxation constant in the fully corrected expressions originate from the 
summation of infinite number of diagrams; in the well-known result of Ohmic BO model 
the bilinear coupling between the system-bath is fully taken into account; this is the exact 
prediction from a simple reasonable model and we concern with the relaxation of fully- 
dressed states in the exact result of BO model. On the contrary, the relaxation constant in 
Eq. flT?D, is the result of the same model but with the second-order (in the coupling strength) 
approximation. Nonetheless in some context the second-order result has been favored while 
the full-order result has been questioned. |33|, [45 



As we show below we can distinguish the above two models ((Jl6|) or (|T7])) by some two- 
dimensional experiment by checking existence or absence of certain peaks. In other words, 
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whether the coherence (off-diagonal) relaxation constant which depends only on the quantum 
number difference (where r m+n m = r n ) and the level independent population relaxation 
is appropriate (as the first-order picture) or not might be checked experimentally. 

Note that if the system has some sort of anharmonicity such as the anharmonicity of po- 
tential [ 1A ] or the nonlinear system-bath coupling , the relaxation constants do not hold 
the relation I^i = Tio etc., even we take into account higher-order system bath interactions. 
Then the number of Liouville paths involved in the optical processes increase dramatically 
especially when the system-bath interaction is very strong. Also if the laser-molecular in- 
teraction is much shorter than the time duration of the system-bath interactions, one has 
to regard the relaxation rate as a function of time, i.e. T nm (t). In such case, the equation 
of motion approach is more appropriate than the diagrammatic approaches, although it 



requires computationally expensive calculations. [47, 48, 49. 50 [46 [51 



We comment on confusion in the literature with regard to the Redfield theory, one ex- 



ample of which is Eq. (17). The Redfield theory without the rotational wave approximation 
(RWA) is equivalent to the Fokker-Planck equation. [[46 



The time evo- 
lution operator in the Liouville space from the state \k,l)) = \k) (/[ to = (i\ ■ ■ ■ \j) 
is then expressed as ((i,j |e -1 ^ /)), where H x is the quantum Liouvillian and f 
is the damping operator (Redfield operator). In energy-level representation, \k,l)) is the 
eigen-function of the Hamiltonian but not the eigen function of T, which makes difficult 
to evaluate this propagator. However, one sometimes "reads off" the damping constant 
directly from the Redfield tensor elements T^f-i and incorporate them in the propagator 
as ((i, j | e^ 1 ^ 1 ~ r ^ kl ) t | k, /)), which can not be justified from the coordinate representation 
model. |33|, |44| Accordingly, this ad hoc methodology possesses a flaw in the sense that the 
theory thus obtained does not converges to analytical perturbative results such as obtained 
by the Brownian oscillator model. It is possible to evaluate effective tensor element rjjjjy by 
solving the equation of motion such as the Fokker-Planck equation with linear and nonlinear 



system-bath interactions, [46, 47, Hq, 49, 50, 51] but the calculated results are quite different 



from the Redfield tensor elements. |46| 
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VIII. MULTI-MODE SYSTEM 



Extension to the multi-mode system, whose characteristic modes are represented by {Q s }, 
{M s }, and {7s}, is straightforward. |2^, 0|]l7|, |35[ We expand the dipole or polarizability 



operator as 



and we denote the Liouville state by 

IK}) (Ml = (K> Kl)i ■ • ■ (K) (K\) s • • - , 

where {n s } = (ni,n 2 , ■ ■ ■) is the quantum-number of the corresponding mode. Here and 
hereafter, we use the notation in which \n s , n s i) (m s ,m s i\ stands for the state where the 
mode s and s' are in the states \n s ) {m s \ and \n s >) (m s >\, respectively. For example, |0, 1) (2, 3| 
means that the first and the second modes are in the ground and the first excited ket states 
while they are in the second and the third excited bra state, respectively. 

The factor (in the Feynman rule) for the transition is well explained by example. The 
transition, 

|0,0) (0,01-12,1X0,01 

caused by the operator (Q^) 2 is associated with the factor 

l(4 112) +4 m) +4 211) )(2,l|(Q (1) ) 2 Q (2) |0,0)/3! = g 112 '^^! while 
the transition (again caused by (Q (1) ) Q (2) ), 

|0,0) (0,0| — |0,0) (2,1| 

is associated with the same factor with the minus sign. If the above transition occurs at the 
last time, however, we have to omit the factor i/h as in the single-mode case. 
Note here that the transition of the type, 

|0,0) (0,0|-|1,0) (1,0| 

cannot occur at once, but 

|0,0) (0,0| - |1,1)(0,0| 

can occur; bra and ket excitation can never occur simultaneously, that is, the simultaneous 
multi-transition can occur exclusively for the ket state or for the bra state. 
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The time propagation factor of each mode in the state (\n) (m\) s during a (positive) time 
duration t is given by Q- i \ ri ~ m ) n ' t for the harmonic system without dissipation. 

In the multi-mode case, the diagram explicitly written in the square bracket D2 in Fig. 
P represents either a single-mode process, 

|1) (0| -> |1) (1| (mode s) 
► - (mode s') 

where — implies no time propagation, or a two-mode process, 

|1) <0| - |1) <0| (modes) 

> ( i8 J 

> |0> <1| (modes') 

which is explicitly shown in the square braket D2 in Fig. [7[ In other words, in the multi- 
mode case, quartet D2 in Fig. || represents the quartets displayed in Fig. ||. 

By use of the above rules in the multi-mode case, we see that the propagator of the 
process in Eq. ([TJJ) is given by e~ lQsTl • e - l ( n "-n s ,)T 2 Decause 1 1 ; 0) (0, 0| propagates for 
T\ and |1,0)(0, 1| for T 2 . The remaining interaction factors s) (1,0|QW|0,0) • 

(-|) x[ s>) (0, 0| g( s ') |0, 1)- (x { * s>) + x ( / s) ^ (0, 1| Q (s) Q (s ' } |1, 0) /2! (and the factor 1/2 to avoid 
double counting). Taking into account the other elements of the quartets, we obtain the 
total contribution D2 of Fig. ^| in the multi-mode case in a form: 

,2 

D2 = 4 Re 



ic„,(i) xf'.^ (19) 



^ c -jn a Ti e -i(n s -n s , )t 2 



2.1 /J). 2M S /Q S / 



where c ss / is 1 and 3/2 for s ^ s' and for s = s', respectively. Comparing this with diagrams 
we learn that we should associate \n s ,n s >) (m s ,m s i\ with Qn]m s + flnjm a ,- These 4 quartets 
correspond to 4 diagrams in Fig. |] (in the dissipation- less case). 
In this way (taking into account the effect of dissipation), we have 

i2 (2) (T!,T 2 ) (20) 
= & + H * + B s + C s + Dl s + D2 2 ) 

s=l,2 

+ ^ A2 SS / + ^2 (B ss ' + C ss / + D ss /) , 
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FIG. 7: Two-mode processes. It should be noted that there are no counterparts of I', II', Al, Dl. 
where the prime in the expression, ^2 S , , implies that the terms with s — s' are excluded 
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FIG. 8: Quartets represented by the quartet (D) of Fig. 
in the sum. Here, each term is given by: 

I, = -/ M e- r » ri ~ r io T » cos {2Q S T X + ( S T 2 ) 
Us = f„ e - r »^- r ^ cos (2C s Ti + CT 2 ) 
A2 SS , = -f ss , e ~ T ^- r ^ cos (C S T X + Q S ,T 2 ) 
B^A.e-^-^oosCC.rx-C.Ta) 
B ss> = / ss , e - r io Tl - r oi° T2 cos (CTx - ( S ,T 2 ) 
C s = -f ss e-^-^ cos (CTj + 2C S T 2 ) 

<^ss' — —Jss' e V 7 

xcos (c.Ti + (& - GO r 2 ) 

Dl s = -i/ ss e- r ^-^ 2cos(CsTi) 
D2. = |/„e- r ^-^cos(C^i) 

D2 SS / = -/ M /e V 7 

xcos (^ + (0-0)^) 



with 



fs 



(s) (s') {as') 
^ 1 ^ 1 *^ 2 



2M S ( S M S ,Q 

We remark the following: (1) I ss / and II SS / always cancel out while l s and II S cancel out only 
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if 



r( s ) — 

1 10 — L 21 ■ 



(2) The sum (Al + A2) s is just given by setting s' — > s in A2 SS /. (3) When we put 

rW= 7a /2for |n s ) (n s \ 

r mn = K ~ m s \ 7 s /2 for \n s ) (m s \ (n s ^m s ), (21) 

the above expression reduces to the result of the fully corrected Brownian oscillator model. 
If we employ the model with the relaxation constant, 

(s) = n + m 

mn 2 ' s 

this leads a different result; one of the feature is the survival of the single mode terms I s and 
II,. 

IX. FEYNMAN RULE IN FREQUENCY DOMAIN 

In the frequency domain, we study the quantity 



poo roc 

/ dcui / duj 2 e^ Tl+ ^R^{T x ,T 2 ) 
Jo Jo 



The frequency domain expression is obtained by using the above propagators in frequency 
domain (or, instead, directly by Fourier transformation of Eq. (EH)). The general propa- 
gating factor in the multi-mode case, e - rT i-* OT i . e -r'T 2 -iQ'T 2 ^ j n ^ e f re q UenC y domain, 

replaced by 

i i . . 
• (22) 

X. 2D SIGNAL FROM EACH LIOUVILLE-SPACE QUARTET 

In this section, we present two-dimensional signals from each Liouville-space quartet 
separately in the fully corrected Brownian oscillator model. In the frequency domain, since 
the signal is a complex number, we show the absolute value of the signal. In the time 
domain, the signal is real, which is directly shown. 
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FIG. 9: Contour plot of the signal from the system with a single mode with weak damping. The 
upper four plots correspond to the separate contribution from each Liouville-space quartets. The 
bottom plot is the sum of them, i.e., the total signal. 

A. Frequency Domain 

1. Single weakly-damped mode 

Fig. ^] shows signals from the system with a single mode (Q = 1, 7 = 0.1, in arbitrary 
unit). Signals from each Liouville space quartet are separately shown. We can interpret 
each peak in the following way: the process represented by \n) (m\ — > \n') (m'\ imply that 
the system is in the state |n) (m| for 7\ and \n') {m'\ for T 2 ; we assign Q nm and T nm for 7\ 
and Q n 'm' and T n i m i for T 2 . This can be symbolically written as 

(Trim; Tn'm') 
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|n) (m\ — > \n') (m'\ 



Actually, the process \n) (m\ — > \n') (m'\ corresponds to the peak at the position at (u>i, 002) = 
(Q nm , Qn'm') with the width in the c^i-axis and u^-axis given by T nm and T n i m i, respectively. 
This results from the expression in Eq. ( p2|) and can be confirmed numerically as we see 
below. 

We note here that we need not consider the contribution from the quartets I and II 
because they cancels out with each other in the fully corrected Brownian oscillator model. 

Quartets A=A1+A2: this be symbolized by |1) (0| — > |1) (0| and its complex conjugate 
|0) (1| — > |0) (1|. The former process can be symbolically written as 



|1><0|-|1> 



(r 10 ,r 10 ) ^ \ (7/2,7/2) ' 

This suggests a diagonal peak (001,002) = (£l,Q) whose widths in the a>i-direction and 002- 
direction are both 7/2; this peak shows symmetric pattern with respect to the two axis, 
which can be seen in the contour plot in Fig. [3[ With the complex conjugate process 
|0) (1| — > |0) (1|, we associate 



|1) <0|-|1> (0| 



(n 01 ,n 01 ) = I -(n,n) 
(r 01 ,r 01 ) ^ \ (7/2,7/2) ' 

Namely, the quartet pair A corresponds to two symmetric diagonal peaks at (001,002) 
±(Q, tt) (see the top left plot of Fig. g). 

Quartet B: symbolically, the association is as follows: 



|1> (0|-|1) (2| 
and its complex conjugate 

|0> (1|-|2)(1| 



(Q 10 ,Q 12 ) = I 

(r 10 ,r 12 ) ^ (7/2,7/2) 
(r i,r 21 ) ^ I (7/2,7/2) 



Namely, we have two symmetric diagonal peaks at (001,002) = ±(fi, — O) (see the top right 
of Fig. @). 

Quartet C: in the similar way, from the association 



|1> <0|-|2) (0| 



(n 10 ,n 20 ) I (n,-2fi) 
(r 10 ,r 20 ) ^ \ (7/2,7) 
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and its conjugate, we should have two significant overtone peaks at (c^i,^) = ±(0,20) 
whose width in the cui-direction is one half of that in the ^-direction; the peak is elongated 
in the second axis as can be seen in the contour plot in Fig. |9| (see the middle left of Fig. 



Quartet D=D1+D2: from the association, 



D2: |1)<0|-|1)<1| 





Dl: |1><0|-|0)<0| 



and their complex conjugate, we should have two significant elongated axial peaks at 
(ux,U2) = (±0,0) (see the middle right of Fig. H). 

The total signal displayed at the bottom of the Fig. ^ shows 8 significant peaks; now that 
we completely know from which Liouville-space path each peak originates, we can assign 
each peak with distinct Liouville-space paths by the following table. 



quartet 


peak positions 


in (uj\ , 0J2) plane 


(A) 


(0,0), 


-0,-0) 


(B) 


(0, -0) 


,(-0,0) 


(C) 


(0,20) , 


-0,-20) 


(D) 


(0,0), 


(-0,0) 



In Fig. H we notice that peaks from quartets (C) and (D) are elongated in the second 
axis. This point is also understood in the above argument, from which we have the following 
table. 



quartet 


width of peaks for (0)1,002) 


(A),(B) 


(7,7) 


(C),(D) 


(7,27) 



2. Double modes (weak damping) 



Fig. [10] shows signals from the system with two weak damping modes (Oj = 1, 71 = 

1, i.e., 



O.lOi, 2 = 0.5, 72 = 0.1O 2 , in arbitrary unit, with the assumption, x\ 
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FIG. 10: Contour plot of the signal from the system with two weakly-damped modes. 

fss' = (CsCs') 1 -)- Signals from each Liouville-space quartet are separately shown. We can 
interpret each signal in the following way. 

Top-left plot of Fig. [H| Two-mode quartet A2 in Fig. [7] is associated with 



|1)<0|--- (modes) _ | (oS,^? 

— > li> <o| (modes') ^ \ (r$,r 



10 > 1 10 

and its complex conjugate; this quartet produces the four cross peaks at (wi,^) = 
=b(Qi, 2 ), ±(J7 2) ^i)- The remaining four diagonal peaks at {u^Uq) = Ox) an d 

±(^2,^2) originate from the single-mode quartets A2 and Al in Fig. |6|, which corresponds 
to the process 

|1) (0| -> |1) (0| (mode s) 
> - (mode s') 

and its conjugate. The widths in the cui-direction and ^-direction for the peak at uj 2 ) = 
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are T^q and r[g\ respectively. In the fully-corrected BO model they are 7s/2 
and 7s'/2 ; respectively. Although there exists the effect of interferences, the relative size 
of the width is consistent with this indication. For example, this is the reason the peak at 
(1,0.5) and (0.5, 1) are elongated in the co>i and a> 2 axes, respectively. In summary, in the 
fully-corrected BO model, the positions of peaks and two-component of widths are given by 



A2/A1 (single-mode): 
A2 (two-mode): 



±(fil,fil) With (7l/2,7l/2) 

±(fl a ,fi a ) with ( 72 /2, 72 /2) 

±(0!,Q 2 ) With (7l/2,72/2) 

±(fi a ,fii) with (72/2,7i/2) 
Top-right plot of Fig. [K| Single-mode quartet B in Fig. | and two-mode quartet B in [7] 
are associated with 

|1) (0| - |1) (2| (mode s) ^ ( (^S^i?) 

— — (mode,') ^\ (r« rg?) 



[1> <0I — > — (modes) I {^it^oP 




— - |o> <i| (mode,') { (rg,rS°) 

The single- mode quartet produces the four diagonal peaks in the top-right plot, while the 
two-mode quartet the four cross peaks. The widths in the two directions for the diagonal 
peaks are given by (t^q , T^ 2 j while those for the cross peaks by (t^q, • in summary, 



we have 

B (single-mode) 



B (two-mode): 



±(fii,-fli) with (7i/2,7i/2) 

±(n 2 ,-n 2 ) with (72/2,72/2) 
±(fii,-n 2 ) with (71/2,72/2) 

±(fi 2 ,-^i) with (72/2,71/2) 

Middle-left plot of Fig. [U]: Single-mode quartet C in Fig. | and two-mode quartet C in 
[7] are associated with 

|1> <0| - |2) <0| (modes) f^SM?) 
— (modes') ^\ (rg^ff) 

|1) (0| - |1) (0| (mode s) f (n$, + ft^) 

- - |i) (o| (mode s') ^ \ (rg, r^> + 4 s ? 
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C (single-mode): 
C (two-mode): 



The single-mode quartet produces the four overtone peaks in the middle-left plot while the 
two-mode quartet the four cross peaks. In summary, we have 

±(n 1 ,2Q 1 ) with (71/2,71) 
±(ft 2 ,2ft 2 ) with (72/2,72) 

±{Q U Q 1 + Q 2 ) with ( 7l /2, ( 7l + 72 ) /2) 
±(Q 2 , Q 2 + fii) with ( 72 /2, ( 7l + 72 ) /2) 
Middle-right plot of Fig. [K]: Single-mode quartets Dl and D2 in Fig. |6] are associated with 

|1> (0| - |0) (0| (mode s) ( (fift?,^) 

— — (modes') ^\ (r« r«) 

|1) (0| -> |1) (1| (mode s) f (^S^i?) 

— — (modes') ^\ (rw.rj?) 

while two-mode quartet D2 in |7| are associated with 

|1> (0| - |1) (0| (mode s) f - fff) 

— Ho)(i| (modes') ^\ (rg, rg + r^) 

The single-mode quartet produces the four axial peaks in the middle-right plot while the 
two-mode quartet the four cross peaks. In summary, we have 

±(fii,0) with (71/2,71) 
±(fi 2 ,0) with (72/2,72) 



D1/D2 (single- mode): 



D2 (two- mode): 



±(fii,fil-fi 2 ) with (f ,2l±2i; 

±(O 2 ,O 2 -Q0 with (f ,2i±a; 



Note here in the fully-corrected BO model, we have Tq = T\i so that the widths from the 
single-mode quartets Dl and D2 are the same in the above. 

The total signal is displayed at the bottom of the figure; we can assign each peak with 
distinct Liouville-space paths or energy-level diagrams (as in Fig. |8|). 



B. Time domain 



Figs. [XT] shows the contour plots of peaks from each quartet for a single over-damped 
mode system. Each quartet contributes to the total signal rather different way. This suggests 
the possibility of Liouville-space-path selective spectroscopy. 
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FIG. 11: Contour plot of the signal from the system with a single over-damped mode. 

XI. SIGNALS FROM BROWNIAN OSCILLATOR MODEL AND REDFIELD- 
TYPE MODEL 

In Fig. 12, we compare results from two models: (1) Brownian oscillator (BO) model (the 
system-bath interaction is fully taken into account) where we put Eq. ([16]). (2) Redfield-type 
model (RT) where we put Eq. (|H]) with the replacement ( s — > Q s (no frequency shift). 

Top: the right plot from RT model has extra peaks at on the left (BO) at (u^,^) = 
±(2, 1). They originate from the survival of the quartets I and II in Fig. |6|. 

Middle: on the left plot (BO) there exist extra peaks at (cji,^) = ±(1 ; ~ !)• This 
corresponds to the single-mode quartet B in Fig. || For this process, the relaxation constants 
associated with the u;2-axis, T^, in BO and RT are given by 7 and 37, respectively; the 
relaxation in RT is much faster, which explains the disappearance of the peaks. The peaks 
at (a>i, ^2) = ±(0.5, —0.5) still survives because these peaks not only come from the single- 
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Brownian Oscillator m odel 



Redfield-type model 




FIG. 12: 2D signal from the two models. Top: single mode, f2i = 1. Middle: two modes, f2i = 1, 
= 0.5. Bottom: two modes, Oi = 1.2, 0,2 = 0.3. Depending on parameters, the difference 
between the models manifests as existence or non-existence of diagonal peaks. 

mode process B: in this case, the peaks from quartet I and II overlap with those from other 
quartets. 

Bottom: on the right plot (RT) exists extra peaks at (cui,^) = ±(0.3,0.6). They corre- 
sponds to the survival of I and II in Fig. || 

In summary, the detailed situation depends on parameters. However, they have one thing 
in common; the difference between the models manifests as existence or absence of certain 
peaks. In the numerical results given above are all in the weak damping regime (7 ~ 0.1Q). 
The weak effect, nonetheless, affects the existence and absence of certain peaks. This is 
because the damping constants directly matter in the cancellation mechanism of certain 
processes. Note that the situation is completely different for weak potential anharmonicity 
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or nonlinear polarizability. Such weak effects, on the contrary, do not concern delicate 
cancellation mechanisms. 

If the system exhibits non-weak anharmonicity of potential or the nonlinear system-bath 
coupling, as mentioned before, there may be the peaks at the similar position predicted by 
the Redfield-type model. Such mechanism, however, affect not only the existence of these 
peak but also the entire profile of signal, which involves different Liouville paths. The careful 
study of the signal in frequency domain shall be the critical test of the Redfield-type model. 



XII. CONCLUDING REMARKS 



We stated an interpretation of the energy-level diagrams in the Liouville space and sum- 
marize the relationships between several diagrammatic representations. We emphasized all 
the diagrammatic representation reduces to unique interpretations in Liouville space, via 
which we can write down analytical expression by a Feynman rule. 

We have given examples in which each Liouville process make distinctly unique contri- 
bution to two-dimensional signal; the selective detection of quantum process by ultrafast 



spectroscopy might be possible, for example, by utilizing the phase matching condition. |37 



By suitably prepared spectroscopic configuration, we might be able to concentrate on a cer- 
tain quantum processes, which allows simpler analysis and more quantitative understanding. 
Such Liouville-space-path selective spectroscopy might be promising. As photon echo can 
be distinguished from the pump-probe via phase matching condition, we could differentiate 
spectroscopic methods by the peaks they produce. 

Energy-level diagram is useful in interpreting the physical process but it is so only after 
confirming the diagram certainly makes a non-zero contribution possibly by other method. 
For example, in the (fully-corrected) Brownian oscillator model, we can assume the ini- 
tial state of the diagram is the ground state; this is because we know that other initial 
states result in the same contribution from a separate calculation. Another example is the 
cancellation of I by II of Fig. Q 

In this respect, diagram in the field-theoretical context, for example, introduced in [[27] has 
some advantage. Number of the diagram to be consider is considerably smaller and analytical 
expression is much simply obtained; in the case of R( 2 \T 1 ,T 2 ), we have only to consider just 
two diagrams in total each being given by the product of two certain propagators. This 
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is because the cancellation is always automatically taken into account in this method and, 
in addition, quartets are summed up from the beginning in a simpler form. However, this 
conceals physical processes in the Liouville space. 

Note here that we have to carefully check out all possible cancellations even if we incor- 
porate the phase matching conditions into the response function, as have been utilized in 
the electronically resonant experiments such as photon echo and pump-probe to detect the 
different contributions of the Liouville paths by choosing the laser wavevectors and frequen- 
cies. [ l] This is because in vibrational spectroscopy such as IR echo, the time durations of 



laser pulses are much shorter than the time periods of molecular vibrations. [[37]] In addition, 
if the initial temperature of the system is higher than the excitation energy of vibrational 
levels (as in the case of low frequency modes), or if the nonlinearity of the dipole or Ra- 
man transitions are important, ||, [l(], |27j we have to include a number of Liouville paths 
especially in higher order spectroscopy; the assignment of the peaks to some Liouville paths 
become nontrivial. 

As for the mechanism of relaxation, we have only considered the system bilinearly coupled 
with bath. We constructed the Feynman rule by starting from the rule in the case without 
damping and then by replacing the propagator so that it causes damping with an appropriate 
choice of the relaxation parameters T mn . One may think that the set {T mn } is an arbitrary 
set of parameters to fit experimental data; in the case of vibrational spectroscopy, however, 
r m „'s have to satisfy certain universal relationships, for example, to satisfy the detailed 
balance condition. In addition, the validity of the rotating wave approximation (RWA) and 
the Markovian approximation associated with the second order perturbation of the system- 
bath interaction might become questionable in vibrational spectroscopy; the characterization 
of the relaxation processes by simple rate constants such as T\ and T2 might not work. Note 
here that, although there are some restrictions, one can calculate the signals without using 



such approximations for Brownian model even in the anharmonic case, [gg, |27], [31], [48| In 
order to verify the consistency of theory, it is important to compare the results from energy 
level models and the Brownian motion model where the latter is based upon a microscopic 
picture. 

In order to demonstrate how approximations for relaxation processes can change the 
results, we presented the 2D signals from the Redfield-type model and (full-order) Brownian 
Oscillator model, and we observed that two models give peaks at different positions even 
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for weak damping. This, in turn, suggests high sensitivity to the damping mechanism of 
2D spectroscopy. This situation is in good contrast with cross peaks associated with mode- 
coupling of anharmonic or nonlinear origin. They might be fairly strong to be observable 
with stronger diagonal peaks. On the contrary, the cancellation mechanism is subtle and, 
thus, weak damping effect can cause a drastic difference. 

One of the purposes of our paper is to bridge the two complementary approaches of the 
coorinate-based and the energy-level based models. The results allow us a useful interpreta- 
tion of the coordinate-based model in the energy-level language. We should note, however, 
that this interpretation becomes precise only in the weak damping limit. Nonetheless, we 
believe that it is useful to have a common interpretation for the two approaches in certain 
situations. 
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APPENDIX: EXPRESSION FOR R^(u u u 2 ) 
2R ( - 2 \ijJi,uj 2 ) is given by Eq. fl2"0p with each term expressed as 
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